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We present estimates of parameters which model the primordial density fluctuations using 
Q ■ CMB data from the VSA, 1st year WMAP, CBI and ACBAR along with external priors from 

2dF and the HST. Taken at face value there is significant evidence for a running spectral index 
when only scalar fluctuations are considered. The significance of this result is increased by 
the inclusion of tensors and, although the result is weakened by the inclusion of 'running of 
the running', the basic trend toward significantly negative running remains. We discuss the 
implications of this result and how it might be strengthened by future observations. 
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1 Introduction 

The Very Small Array (VSA) is a 14 element interferometer which has been operating in the Ka 
band [y ~ 26 — 36GHz) at the Teide Observatory in Tenerife since 2000. The present discussion 
concerns the most recently presented measurements'*™ which are discussed elsewhere in these 
proceedings^. They represent ~ 7000 hours of actual integration on a total of 33 fields by the 
VSA in its extended configuration. The specific results presented here focus primarily on the 
constraints obtained on the primordial power spectrum in conjunction with the 1st year WMAP 
data^, the most recent CBI results^ and those from ACBAR"^. They follow up and expand 
upon those already reported^. 

We use the COSMOMC package to perform our calculations and maginalize over the other 
cosmological parameters Q.\ ) h? , £l m h 2 , h, As and z re as described in previous work™. We use the 
same flat priors except for the redshift of reionization, z Te , for which we use the range (4, 50). 

The scalar component of the primordial power spectrum is modelled as a power series in 
log{k/k c ), 

log P s (k) = log A s + (ns - 1) log(fc/A0 + ln R [log(k/k c )} 2 + ^n RR [log(A;/fc c )] 3 + .. (1) 



Table 1: Estimates of ns, «-r, hrr, nx and At /As for models A, B and C. We have used CMB data from WV 

and WACV with external priors from 2dF and HST. 





CMB 


EXT 


n s 


reR 


™RR 


?1T 


At/As 


A 


WV 
WACV 


HST 
2dF 
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2dF 


0.93 ±0.06 
0.93 ±0.05 
0.91 ±0.06 
0.90 ±0.05 


-0.08 ± 0.04 
-0.05 ± 0.03 
-0.08 ± 0.04 
-0.06 ± 0.03 
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WACV 


HST 
2dF 


1.26 ±0.12 
1.18 ±0.11 


-0.11 ±0.04 
-0.08 ± 0.04 






-0.02 ±0.49 
-0.08 ±0.42 


0.49 ±0.39 
0.27 ±0.22 
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HST 
2dF 


0.89 ±0.07 
0.89 ±0.08 


-0.18 ±0.10 
-0.09 ± 0.09 


-0.04 ±0.05 
-0.01 ±0.05 











where k c is the pivot scale, ^4s is the amplitude of the power spectrum, ns is the spectral index, 
^R ("-run in ref.[2]) is the running of the spectral index and n^R, is the running of the running of 
the spectral index. For this to be considered a valid power series expansion, one will require each 
of the terms to be smaller than the previous one over the range of k relevant to cosmological 
observations (log(A:/A; c ) sa 4). This requires that |ns — 1| 3> 2nR 3> 8nRR/3. 

We will consider three different models (A) Scalar fluctuations only with n^R = 0; (B) 
Scalar fluctuations with n^R = and also tensor fluctuations with power spectrum Pt(^) = 
Ar(fc//c c ) nT ; (C) Scalar fluctuations only, but with n^R ^ 0. For models A and C we will 
use k c = 0.05 Mpc~ and for B, k c = 0.002 Mpc~ . Since we find that models with significant 
running are favoured, the preferred values of ns will appear to be very different for model B 
when compared to A and C. 

In addition to data from CMB experiments, we will use information from external sources 
to augment our analysis. In particular, we will use information from the HST Key Projection 
h = 0.72 ± 0.08 with the errorbars assumed to be Gaussian (denoted HST) and from the 2dF 
Galaxy Redshift survey"^ which constrains the matter power spectrum (denoted 2dF) 
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Figure 1: Marginalized 2D distributions for (a,left) ns — hr and (b,right) wr — z vc for model A using WACV and 
HST. The contours represent the 68% (inner) and 95% (outer) confidence levels. 



2 Results & discussion 



The results of our analysis are presented in Table ^ f° r models A,B and C for the CMB data 
from WV (WMAP±VSA) and WACV (WMAP±ACBAR±CBI±VSA), and the external priors. 
The most striking feature in the case of model A is that wr < is preferred and for the most 



part at > 2a. We see that when one uses the HST prior, the preferred value of riR is lower than 
for the 2dF prior. We note that the preferred values of wr are at odds with the conditions for 
the power series expansion to the valid, but not at the la level. Models with br < have an 
absence of power on both large and small scales relative to scale invariance, with ns > 1 on very 
large scales and ng < 1 on very small scales. It appears that the data favour such models. 

It is interesting to consider what slow-roll inflation predicts for h-r. For a power-law model 
with V{4>) c* 4^i a simple approximation yields ns(N) « 1 — ((p + 2)/2N) (< 1 for all N) 
where N is the number of e-foldings from the end of inflation when the corresponding comoving 
wavelength leaves the horizon. Hence, one can deduce that tir ~ — 2(ns — l) 2 /(p + 2) which the 
results clearly violate. In general all slow roll inflation models predict that tir ~ — (r&s — l) 2 and 
our results present a significant challenge to them. Multi-field inflation models are much less 
well constrained by these results. 

One bizarre feature of our results is that the preferred values of z re are greater than one 
might expect, even in the most baroque models for reionization. For model A and WACV the 
constraint on z Te = 23±7 for the HST prior and z re = 18±6 for the 2dF prior. We have presented 
the degeneracies in the riR — z Te and ns — nR planes for the HST prior in Fig.Q We see that large 
negative values of h.r correspond to large value of z re . A unique feature of the pres ent analysis is 
that we include the full range of values of z re whereas various other analyses^ ' ' use restricted 
ranges. While excluding high values of z re might be suggested by arguments about the nature of 
reionization, leaving them out of the analysis with a flat prior can bias the confidence limits on 
the other parameters by introducing a model dependent slicing of the likelihood surface. This 
issue will be resolved by upcoming WMAP polarization data which should constrain z re much 
more tightly. 

A wide class of inflationary models predict a significant tensor amplitude. Our results for 
model B show that the preference for riR < is more significant when they are included in the 
analysis. This is because models with At/As ^ prefer more negative values of tir than for 
At /As as illustrated in Fig. |^a) for WACV and the HST prior. 

One might expect that if riR is significantly < as suggested by our results then one might 
also expect that hrr ^ also. Our results for model C appear to confirm the preference for 
nR < 0, albeit at less significance, while appearing to be consistent with tirr = 0. There is a 
strong degeneracy between riR and ?t.rr as illustrated in Fig. [2fb). 

We should note our previous work highlighted a number of possible systematic effects which 
could be at work here. The parameter ur can only be constrained by making measurements 
over a range of angular scales making it particularly sensitive to a wide range of effects. In 
particular, we highlighted issues pertaining to the overall calibration uncertainty of the data, 
possible errors in the point source subtraction, and the possibility that the measured WMAP 
quadrapole is contaminated by foreground emission. 

Finally, in order to re-inforce the point that the data favour models with a dearth of power 
on small scales, we have also considered a model with a broken power law spectral index. 
Specifically, we have used n$(k) = n\ for k < k c and ng = n 2 for k > k c , with k c = 0.05 Mpc -1 . 
We find that for WV+HST n t = 1.05 ± 0.05 and n 2 = 0.54 ± 0.22, whereas for WV+2dF 
n\ = 0.98 ± 0.05 and n 2 = 0.87 ± 0.15. Although this should not be taken too seriously as 
physical model, it illustrates the point very strongly. 

3 Future observations using the super-extended array 

It is planned that the VSA will be reconfigured in the late summer of 2004 and fitted with larger 
mirrors to achieve much higher angular resolution. It is expected that it will be possible to 
make observations which are sensitive to £ = 500 — 2500. Given that the primary CMB signal 
(in temperature) at £ = 2000 is a factor of about 3 weaker than at £ = 1000, it is planned to 
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Figure 2: Marginalized 2D distributions for (a, left) na — At /As for model B using WACV and HST, and (b, right) 
Tin — «rr for model C using the same data. Confidence levels as in fig. 

upgrade the array electronics to give a factor ~ 3 improvement in sensitivity to maintain the 
same signal-to- noise. This will be achieved by an increase in the bandwidth from A/ = 1.5 GHz 
to 6 GHz and a reduction in T sys from 35K to 25K. 

We have produced simulated errorbars for an 18 month observing program and used them 
to create a mock COSMOMC file based on a fiducial model which is the best fit to WV+HST. 
We found that the derived error on the running of the spectral index Ah-r reduces by about 30% 
from AriR = 0.038 to AriR = 0.027 suggesting that observations at high £ can have a significant 
impact on the constraining the primordial power spectrum. 
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